How does confinement affect the dynamics of viscous vesicles and red blood cells? 
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Despite its significance in microfluidics, the effect of confinement on the transition from the 
tank-treading (steady motion) to the tumbling (unsteady motion) dynamical state of deformable 
micro-particles has not been studied in detail. In this paper, we investigate the dynamics of a 
single viscous vesicle under confining shear as a general model system for red blood cells, capsules, 
or viscous droplets. The transition from tank-treading to tumbling motion can be triggered by the 
ratio between internal and external fluid viscosities. Here, we show that the transition can be induced 
solely by reducing the confinement, keeping the viscosity contrast constant. The observed dynamics 
results from the variation of the relative importance of viscous-, pressure-, and lubrication-induced 
torques exerted upon the vesicle. Our findings are of interest for designing future experiments 
or microfluidic devices: the possibility to trigger the tumbling-to-tank-treading transition either by 
geometry or viscosity contrast alone opens attractive possibilities for microrheological measurements 
as well as the detection and diagnosis of diseased red blood cells in confined flow. 
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I. INTRODUCTION 

The dynamics of deformable micro-particles (e.g., 
droplets, polymeric capsules, phospholipid vesicles, or 
blood cells) under flow is a fascinating fundamental prob- 
lem with increasing relevance for technological applica- 
tions, for example, in food processing, drug delivery, or 
designing lab-on-chip devices. The present paper con- 
stitutes a contribution to the study of vesicle and red 
blood cell (RBC) dynamics, an attractive subject that 
has led to numerous publications in the last two decades 
(see [l[ and the references therein). A key aim of these 
works is to derive constitutive physical laws to bridge 
the details of the dynamics of each single particle (at the 
microscale) to the global behavior of a suspension (at 
the macroscalc). This bottom- up approach is generally 
used to study the dynamics of so-called complex fluids 
(e.g., blood) since the classical way of treating these as 
a single continuum medium fails to capture the relevant 
physics. It is required especially when the fluid is con- 
fined to length scales comparable to the particle size, for 
instance, blood flow in the microcirculatory system. The 
same scenario is faced in lab-on-chip devices where the 
composed fluid (carrying fluid with suspended particles) 
flows in micro-channels with cross sections of similar size 
as that of the transported particles. Moreover, it was 
found that confinement - due to the presence of the walls 
- alters the dynamics (see @, [|[ for polymers) and stabil- 
ity of deformable micro-particles (see 0t9 for droplets). 
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The present study is devoted to investigate the inter- 
play between the confinement and the viscosity contrast 
(ratio between internal and external fluid viscosities) on 
the dynamics of viscous vesicles and RBCs. This study 
is relevant to elucidate, for example, the challenge of 
the reproducibility and the accuracy of blood viscosity 
measurements that results from the interplay between 
rheometer wall effects and the dynamical state of each 
individual RBC. On the one hand, confinement can lead 
to an overestimation of the measured effective viscosity 
of blood as mentioned in 0. On the other hand, the 
transition of the dynamical states of vesicles and RBCs 
induced by varying the viscosity contrast was found to 
minimize the effective viscosity Q. However, to our 
knowledge, a systematic study that takes into account 
both effects (confinement and viscosity contrast) is lack- 
ing. It is known 0-[l2| that an unconfined vesicle or 
RBC subjected to shear flow undergoes either a steady 
liquid-like motion called tank-treading (the vesicle main 
axis assumes a steady inclination angle with the flow di- 
rection while the membrane undergoes a tank-trcading- 
like motion) or unsteady solid-like motion called tum- 
bling (the vesicle flips as a solid particle [Hj|). These 
dynamical states depend on three dimensionless control 
parameters: (i) the swelling degree A, (ii) the viscos- 
ity contrast A, and (iii) the capillary number Ca (ratio 
between viscous and membrane bending forces) [14l4l7l |. 
The transition between the two dynamical states can be 
induced by varying these three control parameters. In 
this paper, we show that there is an additional control 
parameter: the degree of confinement \- 

In the reported simulations, we place a vesicle in a 
linear shear flow. This way, one avoids the interplay 
between the vesicle deformability and nonlinear velocity 
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profiles which would lead to complex dynamics, even in 
the absence of a viscosity contrast [H, [01 . We keep the 
Reynolds number, the capillary number, the shear rate, 
and the swelling degree unchanged. Consequently, we ex- 
amine how the system reacts to the variation of only two 
control parameters: the viscosity of the enclosed fluid (an 
intrinsic vesicle property) and the degree of confinement 
(ratio of the vesicle size to the height of the channel). The 
effect of confinement on the transition of vesicle dynamics 
was briefly studied by Beaucourt et al. using phase- 
field simulations: For given x, the transition was induced 
by varying A. In contrast, we induce the transition by 
varying only the confinement, keeping all fluid and vesi- 
cle properties (especially A) unchanged. This provides 
the interesting possibility to trigger the transition in ex- 
periments by merely adjusting the wall distance without 
replacing the vesicle or the suspending fluid. Addition- 
ally, we report how the physical quantities characterizing 
each dynamical state vary with confinement. 

II. MODEL AND METHOD 

Vesicles are closed membranes made of phospholipid 
molecules, e.g., dioleoylphosphatidylcholine (DOPC). 
They constitute a good biomimetic model for living cells, 
for example, to study blood cell dynamics. They are also 
used for micro-encapsulation of active materials in drug 
delivery. At room temperature (~ 25°C), most of the 
phospholipid membranes, in particular biological mem- 
branes, are in the liquid state. Their thickness (~ 5 nm) 
is negligibly small compared to the vesicle size (typi- 
cally ~ 10 lira, for giant unilamellar vesicles). Therefore, 
vesicle membranes are considered as two-dimensional in- 
compressible Newtonian liquids. As a consequence, the 
membrane area cannot undergo extension or compression 
since such deformation modes would cost a lot of energy. 
This implies local and global conservation of the vesi- 
cle area, as is the case for RBCs [2(|. Mathematically 
and numerically, this constraint can be achieved by us- 
ing a local Lagrangian multiplier a that plays the role 
of an effective surface tension. The associated tension 
energy is given by Eg = J gQ a(s) ds, where the integra- 
tion is performed over the membrane surface dQ. The 
mechanical deformation mode that costs less energy is 
bending, Eb = ^ J af2 c 2 ( s ) ds, where kb is the bending 
rigidity and c(s) the local membrane curvature. When 
the membrane is bent due to hydrodynamic stresses, and 
therefore brought out of its equilibrium configuration, 
it exerts a reaction force back on its surrounding fluid. 
The force is derived by taking the functional derivative 
F = S(Eb + Es)/8y. In two-dimensions (2D) this results 
in [H| 
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where n and t are, respectively, the unit normal and 
tangent vectors. Our present model can be applied to 



vesicles as well as to RBCs. In 2D, both particles ex- 
hibit similar dynamics [HI, [l9|, [2l|. However, in three- 
dimensional models, we need to account for the RBCs 
shear elasticity. This leads, in some situations, to dynam- 
ics not observed for vesicles (with zero shear elasticity), 
for example, the swinging motion Q. 

The vesicle membrane is a free-moving boundary 
whose shape and position is not known a priori. Its dy- 
namics results from its interaction with the flow of the 
surrounding fluid. In the present work, the membrane 
dynamics is computed using a front-tracking method 
in 2D. The membrane is represented by a moving La- 
grangian mesh (a contour in 2D) and the fluid flow is 
computed on a fixed Eulerian lattice. Instead of solving 
the Navier- Stokes equations inside and outside the vesicle 
directly, we use the lattice-Boltzmann method (LBM). 
The LBM is commonly employed to study th e dy namics 
and rheology of complex fluids (see, e.g., |22h24| ). Both 
the position and velocity spaces are discretized; here we 
use nine discrete velocity directions in 2D. The main 
quantity in the LBM is the density distribution fi(r,t) 
giving the probability to find an elementary portion of 
the fluid at position r with velocity in direction e*. The 
time evolution of fi is governed by the lattice-Boltzmann 
equation with the Bhatnagar-Gross-Krook (BGK) colli- 
sion operator :25| 

fi(r + e t At, t + At)= Mr, t) - — [/,:(r, t) - jf (r, *)] , 

T 

(2) 

where At is the time step, r is a relaxation time re- 
lated to the dynamic viscosity rj via the relation 77 = 
pc 2 (r— 1/2) Ax 2 / At, where p = J2i=o /» IS * ne nm< ^ mass 
density, c s is the speed of sound and Ax the grid spac- 
ing. /? q is an equilibrium distribution. The local fluid 
velocity u = - J2i=o fi e i- Geometrical quantities of the 
membrane (e.g. the local curvature c) required to eval- 
uate the force given by Eq. (fTJ are computed using the 
finite difference method (FDM). The membrane force F 
acting on the fluid is included in the model by adding 
the term (F • e^) to the right- hand-side of Eq. ((2]). The 
two-way coupling between the fluid flow and the mem- 
brane dynamics is accomplished using Peskin's immersed 
boundary method (IBM) (26[. The physical quantities 
computed in each mesh are matched by interpolation. 
More details about the numerical methods can be found 
in Ref. [2lj . In the present work, in contrast to our pre- 
vious article [2l[ , we consider that the viscosities of the 
internal fluid r^nt and of the external fluid ?7 0X t are dif- 
ferent (A = ?7int/??ext 7^ 1)- We locate the fluid nodes 
relative to the vesicle (inside or outside) by the even-odd 
rule [l?) ■ Then the LBM relaxation time is set depending 
on the current fluid node location. Unlike in phase-field 
or level-set methods, the even-odd algorithm does not 
rely on diffuse scalar fields for the viscosity, and no addi- 
tional field equation has to be solved. To our knowledge, 
this is the first time the dynamics of viscous vesicles is 
studied using a combination of the lattice-Boltzmann and 
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immersed boundary methods. 

In the following, a single neutrally buoyant vesicle is 
positioned at mid-distance between two parallel plates. 
A linear shear flow with desired shear rate 7 is generated 
by moving the two plates in opposite directions. All the 
results shown below have been obtained for a Reynolds 
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number of Re = fr/R^/ife 



0.05 (negligible inertia) 



and a capillary number of Ca = ?7cxt7-Ro/ Kb = 0.5 (low 
deformation regime), where R is the effective vesicle ra- 
dius. In 2D, Rq = P/2ir, where P is the vesicle perime- 
ter. In the low deformation regime, the vesicle shape 
is always close to its equilibrium configuration and the 
deformation does not alter the dynamics. The swelling 
degree A = AttA/P 2 (A is the vesicle area) is chosen 
as A = 0.8 in order to compare with the data in [Ioj |. 
which to the best of our knowledge is the only available 
paper where the dynamics of a confined viscous vesicle 
under shear flow is studied. The external viscosity ?7 ex t 
is also fixed in all simulations. In this paper, we vary 
only the two key controlling parameters: (i) the viscos- 
ity contrast A and (i) the degree of confinement x- The 
viscosity contrast A is set to values between 0.5 and 30 
by varying only the internal viscosity rj mi . For healthy 
RBCs under physiological conditions, 7 < A < 13. The 
confinement is defined as x = Rq/W, where W is the 
channel half-height. 



III. RESULTS AND DISCUSSION 

First, we investigate the effect of the viscosity contrast 
A. We compute the physical quantities characterizing 
each vesicle dynamical state when varying A, while keep- 
ing the confinement fixed at x — 0.26. The tank-treading 
motion is characterized by the steady inclination angle 
^ and the membrane tank-treading velocity V, while the 
tumbling motion is characterized by the tumbling fre- 
quency n. Fig. Q] shows the variation of these quantities 
when increasing A from 1 to 16. For convenience, ^ is 
normalized by 7r/6, V and ft are normalized, respectively, 
by the rotational velocity (jRq/2) and frequency (7/471") 
of a rigid cylinder rotating in unbounded shear flow j28| . 
At lower A, the vesicle performs the tank-treading mo- 
tion (left panel). By increasing A, both ^ and V de- 
crease. VP decreases until it vanishes at a critical value 
Ac (here, Ac = 8.25) above which tumbling takes over 
(right panel). In the tumbling regime, Q increases with 
A until saturation at larger A (limit of a rigid particle). 
Even at this finite degree of confinement (x = 0.26), we 
are able to reproduce qualitatively the same dynamical 
behavior known for unconfined viscous vesicles (x = 0), 
reported for example in Refs. |^4l2T|. 

Next, we examine the effect of confinement. We com- 
pute ^ of tank-treading vesicles as a function of A and at 
three different values of x (0.10, 0.26, 0.50). Fig.[2]shows 
the decreasing trend of ^ with increasing A for all three 
degrees of confinement, as seen in Fig. [T] for x = 0.26. 
However, we notice that for both x = 0.10 and 0.26, the 
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FIG. 1: The transition of the vesicle dynamical state 
induced by increasing the viscosity contrast A at fixed 
confinement x = 0.26. Left panel: Both the inclination 
angle ^ and the tank-treading velocity V decrease with 
A. Right panel: The tumbling frequency f2 increases 
with A. Qualitatively, the same dynamical behavior is 
observed for unconfined viscous vesicles (x = 0) (tjjllj. 
All plotted quantities are normalized (see text for 
details). 



inclination angle vanishes at two different values of A. By 
increasing x from 0.10 to 0.26, the threshold of the tank- 
treading to tumbling transition shifts from 6.05 to 8.25. 
For the larger value x = 0.5, ^ decreases as well but 
without reaching zero, therefore, the transition to tum- 
bling motion is not observed. From Fig. [5J we learn that 
confinement increases Ac and even eliminates tumbling 
at higher confinement (for the parameter range explored 
in the present study). 

In Fig. |31 the critical viscosity contrast Ac is plotted 
against the confinement x- This gives the x _ A phase- 
diagram of the dynamical states of a confined viscous 
vesicle subjected to shear flow. It can be seen that the 
transition threshold Ac is pushed up by increasing con- 
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FIG. 2: The inclination angle VP (in radians) versus the 
viscosity contrast A of tank-treading vesicles at three 
different degrees of confinement (x = 0.10, 0.26 and 
0.50). The transition threshold (where "J vanishes) 
shifts to higher values when increasing confinement. 
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FIG. 3: The x-A phase-diagram of the dynamical 
states (tank-treading or tumbling) of confined viscous 
vesicles subjected to shear flow. The confinement shifts 
or even inhibits the transition to tumbling. The range 
of 7 < A < 13 corresponds to healthy RBCs. Typical 
confinement in arteries and arterioles are also shown. 
Circle symbols are data from Ref. [Toj . 



finement. The same trend has been observed by Beau- 
court et al [lO] for 0.1 < x < 0.3. A comparison of the 
present data obtained, by IBM/LBM (square symbols in 
Fig. El), with those reported in [1(| by means of the phase- 
field method (circle symbols in Fig. [3]) reveals a good 
qualitative and quantitative agreement. We consider a 
more extended interval of confinement 0.05 < y < 1.25 
as compared to 0.1 < x < 0-3 reported in [10|. Fig. [3] 
leads to the interesting finding that the transition be- 
tween the tank-treading and the tumbling motion can 
also be induced by varying the confinement x alone, with- 
out changing the viscosity contrast A. For example, by 
taking a tumbling vesicle with A > Ac at lower con- 
finement, we can force it to tank-tread by increasing 
the confinement above a certain threshold xc- We ex- 
pect the same confinement-induced transition (beside the 
shear-induced transition) to happen for RBCs, for which 
7 < A < 13, during their displacement from arteries (low 
confinement: 2W > 100/im, x < 0-1) to arterioles (high 
confinement: 10/im < 2W < lOO^m, 0.1 < x < !)• 

The transition due to confinement occurs via a saddle- 
node bifurcation as shown in Fig. |4l for a vesicle with 
A = 8. The angle is well fitted with a square root law 
~ Vx ~ XC in the vicinity of the bifurcation point xc ■ At 
lower confinement (left panel), the vesicle tumbles. By 
increasing x, the vesicle motion starts to be affected by 
the presence of the walls. The tumbling frequency 51 (tri- 
angle symbols with dotted line) decreases until it vanishes 
at the transition point (xc = 0.249). Beyond this thresh- 
old, and at higher confinement (right panel), the tank- 
treading motion takes over. It has to be stressed that this 
was not expected to happen at x = 0.249. At this degree 
of confinement, the gap between the membrane and the 
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FIG. 4: The transition from tumbling to tank-treading 
of a viscous vesicle (A = 8) induced by confinement x 
(in contrast to Fig. [T] where the transition is induced by 

varying A). Left panel: The tumbling frequency f2 
decreases down to zero with confinement. Right panel: 
The angle ^ increases while the tank-treading velocity 
V decreases with confinement. The dashed-dotted line 
is fit of the angle with a square root law ~ Vx — Xc- 
All plotted quantities are normalized (see text for 
details) . 



wall is three times the vesicle size i?o- This available free 
space would be sufficient for a full tumbling period of the 
vesicle. A physical contact between the vesicle and the 
walls is possible only for x > X*; where x* = 0.75 is the 
confinement for which the vesicle long semi-axis and the 
channel half-height W are equal. Thus, we conclude that 
the tumbling-to-tank-trcading transition induced by con- 
finement is not entirely due to geometrical constraints. 
Fig. f|] (right panel) reveals another astonishing observa- 
tion. The inclination angle "J is found to increase with 
confinement. This result is in contradiction with our pre- 
vious work for non- viscous vesicles (A = 1) where the 
confinement is found to reduce the angle rather than in- 
creasing it [2l|. Furthermore, in Fig. [4] the tank-treading 
velociy V evolves in the opposite way. It decreases with 
confinement. However, this trend is consistent with the 
. The intriguing behavior in Figgis the op- 
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posite variation in ^ and V with confinement Xj unlike in 
Fig. [1] where both quantities decrease with the viscosity 
contrast A. In Ref. plj . the decrease of the inclination 
angle with confinement was explained to be caused by 
a torque exerted by the walls upon the vesicle. Obvi- 
ously, this explanation does not hold for viscous vesicles 
(A > 1). The arising question is then: How does the 
torque incline the vesicle with the flow for A = 1 while 
it does the opposite for A > 1? In both situations, the 
walls, the external fluids, and the membrane properties 
are identical. The only difference is the internal viscosity 
?7i n t and the resulting viscosity contrast A. 

A step forward to answer the above question is to com- 
pare the variation of ^ with x of tank-treading vesicles 
having different viscosity contrasts (A = 0.5, 1, 4, 8). 



5 



9* 

d 

d 

_o 

_d 
U 

d 



71/4 



7C/8 - 



* 

X 


■ 1 ' 

A=0.5 


i 1 
— B — 




A 1 A 


"t> " 




A=4.0 


A 


V 


S^ea A=8.0 


v _ 






w 


i 


i 



0.25 



0.5 0.75 1 
Confinement % 



1.25 1.5 



FIG. 5: The variation of the inclination angle \& (in 
radians) as a function of confinement x of tank-treading 
vesicles having different viscosity contrast (A = 0.5, 1, 
4. 8). For A > 1, \f varies non-monotonically with \- At 
higher confinement, all curves converge to the same 
value "J = 0.23 (radian), x* = 0-75 is the degree of 
confinement when the channel half-height has the same 
size as the vesicle long axis. 



We notice that the vesicles behave differently depending 
on whether A < 1 or A > 1, sec Fig. [SJ For A < 1, 
the angle decreases monotonically with confinement \i 
it decreases rapidly when x > 0.5. For A > 1, the an- 
gle first increases up to a maximum and then decreases 
when the confinement becomes larger. For both cases 
(A < 1 and A > 1), the inclination angles converge to 
the same value (^ = 0.23 radian) at higher confinement 
(x > 1). We observe the same non-monotonic variation 
of the angle with confinement for other values of A > 1 
(not reported here). We are not aware of any work re- 
porting this observation for vesicles or RBCs. However, 
for confined viscous droplets, the increase of "J with x bas 
been observed by Janssen et al numerically [j| and ex- 
perimentally Q , but no physical explanation for this be- 
havior was provided. The authors rather used the angle 
variation with confinement to explain the droplet stabil- 
ity against break-up. Unlike viscous vesicles and RBCs, 
droplets are known to perform only the tank-treading 
motion, and therefore no dynamical state transition is 
observed for them. We explain the observed increase of 
^> with x, for viscous vesicles, by analyzing the pressure 
field. Fig. |5] shows the variation in the pressure around 
a sheared viscous vesicle (A = 8) at different degrees 
of confinement (x = 0.1, 0.5, 1.0). The pressure val- 
ues reported in Fig. [5] are the relative deviation from the 
pressure of the external fluid in the absence of the vesicle. 
Regions with higher pressure are shown in red, those with 
lower pressure in blue. At lower confinement (x = 0.1, 
Fig. UtO) t ne rotational component of the imposed shear 
flow induces a viscous torque in clockwise direction re- 
sponsible for the tumbling motion of the vesicle. Dashed 
closed lines in Fig. [S^ represent the vesicle shapes (taken 
at equal time intervals) to show the tumbling motion. 
The vesicle rotates slowly and spends much time in the 
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FIG. 6: The pressure variation around the same viscous 
vesicle (A = 8) subjected to shear flow at three different 
degrees of confinement: (a) lower x = 0.1, (b) moderate 
X = 0.5, and (c) higher x = 1-0. The red-colored regions 
correspond to higher pressure, while the blue-colored 
ones correspond to lower pressure. In (a) the vesicle 
tumbles. The dashed lines are the vesicle snapshots 
taken at different times. In (b) and (c), the vesicle 
tank-treads. The circle (dotted line) shows the required 
space for tumbling. 



region of — 7r/4 < < 7r/4. By increasing confinement 
to a moderate value (x = 0.5, Fig.[fjb), the pressure vari- 
ations increase. Due to the pressure distribution around 
the vesicle (decreased pressure in the top right and bot- 
tom left quadrants and increased pressure in the top left 
and bottom right quadrants), a counterclockwise torque 
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is induced. This pressure-induced torque competes with 
the torque due to the viscous stress mentioned above. At 
this confinement (% = 0.5), the sum of the two torques 
vanishes. This results in a steady inclination of the vesi- 
cle, instead of tumbling and it explains the occurrence of 
the confinement-induced transition as observed in Figs. [21 
[3] and [4] Increasing the confinement further amplifies the 
pressure-induced torque and consequently increases the 
inclination angle. This explains the increase of the an- 
gle with confinement observed in Fig. [4] and 03 At larger 
degrees of confinement (x = 1.0, Fig. [Bh) the vesicle ap- 
proaches the wall. In this situation, the gap between the 
membrane and the wall becomes thin. However, the vesi- 
cle never touches the walls due to lubrication forces that 
induce a clockwise torque pushing the vesicle away from 
the walls. This coincides with observing a negative local 
wall shear stress (like in [2l| ) . At the largest confinement 
studied (x = 1-25), the steady angle is ^ = 0.23 (radian), 
regardless of the viscosity contrast (see Fig. [SJ . At this 
degree of confinement, the angle is mainly governed by 
the geometry (channel height and vesicle shape) , and the 
other parameters, especially the viscosity contrast, do 
not play a significant role. This is caused by the balance 
of viscous-, pressure- and lubrication-induced torque. 

After explaining the mechanism behind the behavior 
of the vesicle angle with confinement, now we analyze 
the effect of confinement on the way the vesicle tum- 
bles. In addition to the tumbling frequency (shown in 
Fig. [4]) we analyze the phase portrait, where the normal- 
ized angular velocity d^/'ydt is plotted as a function of 
the instantaneous inclination angle ^(t), sec Fig.[7J The 
angular velocity is negative since the vesicle performs a 
clockwise tumbling motion. The value <Ml/dt = —7/2 
corresponds to the angular velocity of a rigid cylinder ro- 




ta) (b) 

FIG. 7: The angular velocity d^/^dt) versus the 
instantaneous inclination angle ^ during the tumbling 
motion of a viscous vesicle (A = 8). (a) the data-set for 

X = 0.1 (square symbols) is fitted perfectly with the 
function A + £?cos(2 , I') (solid line), where ^4=1/2 and 
B = 0.47. (b) the angular velocity versus the angle for 
different degrees of confinement (x = 0.10, 0.15, 0.23). 

At higher x, the angular velocity approaches zero. 



fating in unbounded shear flow [28[ . It has the same value 
independent of 'J. For deflated vesicles (here A = 0.8), 
d\I//dt > —7/2 when the angle is between — 7r/4 and 
7r/4, thus, the vesicle is slower and spends more time in 
that angular region. For the remaining range of angles, 
d\I//dt < —7/2, and the vesicle rotates faster. The solid 
line in Fig. [7^, is a fit of the numerical data-set (square 
symbols) with the function 

1 d* 

— - = A + Bcos(2*) (3) 
7 dt 

introduced by Jcffery to describe the rotation of rigid el- 
lipsoidal particles in unbounded shear flow (l3| . Here, 
A = —1/2 and B is a function of the particle aspect 
ratio (ratio between the long and the short axes in the 
shearing plane). Eq. [3] has since been extended by Keller 
and Skalak Q to take into account the motion of fluid- 
filled particles having a liquid membrane (such as RBCs 
or vesicles). In this case, B depends on the viscosity 
contrast A. The theory of Keller and Skalak (KS) pre- 
dicts the dynamics of a vesicle in unbounded shear flow, 
depending on the value of — A/B. For —A/B > 1, the 
vesicle tumbles, whereas for < —A/B < 1, the vesicle 
tank-treads. It was found experimentally that, for larger 
deformation (Ca ^> 1), the terms A and B deviate from 
their respective expressions given by the theory of KS 
plj . In the present paper, Ca = 0.5 is sufficiently small. 
Therefore, A is considered to be a constant and equal to 
— 1/2. Only B is assumed to be a function of A and Xi 
and used as a fitting parameter. Obviously, Eq. ([3]) accu- 
rately describes the evolution of the inclination angle for 
X = 0.1, see Fig. [7jt. The angular velocity dvP/^di) is 
plotted as a function of ^> for three values of confinement 
(x = 0.10, 0.15, 0.23) in Fig. [7b. For convenience, only 
the region [— ir/8, ir/8] x [—1/4,0] is shown. By increas- 
ing Xi the maximum of the angular velocity approaches 
zero. This means that the motion of the vesicle slows 
down in the region [— 7r/4, 7r/4] where it spends most of 
the time. As a result, the tumbling is delayed. Presum- 
ably, a further increase of x (x > 0.23) would lead to 
a positive value of d 1 ^ /dt and therefore to an unstable 
counterclockwise rotation. Instead, we observe the tran- 
sition to the tank-treading motion. 

By fitting the data in Fig. [TJd with Eq. ^ we find that 
B is an increasing function of confinement x- I n Fig- HI 
we report —A/B versus x, for A = 7, 8, 12. Square sym- 
bols correspond to the values obtained when the vesicle 
tumbles (using a fit as described above). For these points, 
we observe —A/B > 1. The value of — A/B decreases 
until reaching unity at a critical confinement xc where 
the transition from tumbling to tank-treading occurs. 
The corresponding values of xc for each A are in good 
agreement with the phase-diagram in Fig. [3J The trian- 
gles in Fig. [5] denote the values of —A/B computed for 
tank-treading vesicles via the equation —A/B — cos(2$) 
which is Eq. (J3J for a steady inclination angle (d^ /dt = 
0). For these points, < —A/B < 1 holds. As mentioned 
above, in the KS theory, it is sufficient to study the value 
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FIG. 8: Variation of the Keller and Skalak parameter 

—A/B with confinement \ f° r different viscosity 
contrast (A = 7, 8, 12). For each A, the transition from 
tumbling to tank-treading occurs exactly when 
—A/B = 1. The corresponding critical value of x> f° r 
each A, agrees with the one reported in Fig. [31 

of —A/B in order to predict the vesicle dynamical state 
in unbounded shear flow. Interestingly, this also applies 
here to catch the transition, although the system is con- 
fined, while the KS theory does not include confinement 
effects. We emphasize again that here the transition is 
induced by varying \ alone. The viscosity contrast, the 
swelling degree, and the capillary number are all kept 
constant. Thus, the transition from tumbling to tank- 
treading due to confinement can also be explained based 
on the theory of KS. The transition takes place when 
—A/B = 1. Above, we used values of B obtained from 
our simulations. Even though further theoretical study 
is needed to derive a closed analytical expression for the 
X-dcpcndcncc of B, we show that B is an increasing func- 
tion of x- 

IV. CONCLUSIONS 



with A > 1. The confinement-induced transition and as- 
sociated vesicle dynamical behavior result from the com- 
petition between viscous- and pressure-induced torques 
exerted upon the vesicle. At larger degrees of confine- 
ment, the lubrication-induced torque becomes dominant 
and pushes the vesicle to align with the walls. The dy- 
namics of vesicles are well described in the limit of the low 
deformation regime by the theory of Keller and Skalak [9| . 
However, it is limited to unbounded shear flow. Here, we 
predict that the term B in this theory would be a function 
of confinement. The tumbling-to-tank-treading transi- 
tion is also observed for RBCs. It is usually attributed 
to the increase in the shear rate 7, for instance, during 
the displacement of RBCs from large arteries (where 7 
is lower) to the branches (where 7 is higher). A RBC 
also stops tumbling when it is squeezed in a capillary, 
because of the lack of available space. It is necessary for 
RBCs to switch to tank-treading mode in order to pass 
through narrower capillaries. Based on our findings, a 
RBC may be expected to undergo this transition even at 
constant shear rate and at lower confinement (x < 0.5), 
compared to x > 1 encountered in capillaries. The ef- 
fect of confinement can be seen as an extra contribution 
(besides the change in shear rate) to the transition from 
tumbling to tank-treading of RBCs. This transition is 
found to be a contributing factor for shear thinning of 
blood Q. Therefore, since we could induce this dynam- 
ical transition by confinement we expect that this will 
also lead to shear thinning for a suspension of viscous 
vesicles or RBCs. Under these circumstances, the walls 
play the role of inducing the transition; in contrast to the 
Fahracus-Lindqvist effect where the walls induce cross- 
stream migration of RBCs. We expect that our study 
will stimulate further studies, similar to Ref. [29j|, for con- 
fined suspensions of deformable micro-particles in order 
to gain more insight into the link between confinement 
effects, dynamical transition and the effective viscosity. 



Using computer simulations we show that by only 
squeezing a viscous vesicle one can induce the dynami- 
cal transition from solid-like motion (tumbling) to liquid- 
like motion (tank-treading). Moreover, we find that 
confinement slows down the tumbling motion and in- 
creases the inclination angle of tank-treading vesicles 
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